Dynamics of rotating spirals in agitated wet granular matter 
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Abstract. Pattern formation of a thin layer of vertically agitated wet granular matter is investigated experimentally. Due to 
the strong cohesion arising from the capillary bridges formed between adjacent particles, agitated wet granular matter exhibits 
a different scenario compared with cohesionless dry particles. Rotating spirals with three arms, which correspond to the kinks 
between regions with different colliding phases with the vibrating plate, have been found to be the dominating pattern [1]. 
From both top view snapshots and laser profilometry methods, the rotation frequency of the spiral arms is characterized with 
image processing procedures. Both methods reveal that there exists a finite rotation frequency v r at a threshold vibration 
acceleration, above which v r increases linearly with the peak vibration acceleration with a slope strongly dependent on the 
vibration frequency. 
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From natural phenomena such as the formation of 
our galaxy, to the industrial processes such as coating, 
pattern formation is ubiquitous in the universe [2]. Ex- 
plaining pattern formations, especially the accompany- 
ing length and time scales, helps us gain insight into the 
pattern forming systems, which are typically far from 
thermal equilibrium. Driven granular matter, as large ag- 
glomerations of macroscopic particles interacting with 
each other dissipatively, is such a nonequilibrium sys- 
tem [3]. Continuous energy injection is necessary in or- 
der to drive granular particles into a fluidlike state in or- 
der to process or transport them in various industries [4] . 
Even though a all-encompassing theory describing the 
dynamical behavior of granular matter, especially toward 
the limit of dense granular flow, is still far from estab- 
lished, the rich pattern forming nature of granular matter 
provides us the opportunity to gain insights into such a 
nonequilibrium system [5]. 

Besides the cohesionless dry granular matter, more 
and more attentions have been paid to the dynamics of 
wet granular matter recently [6]. Due to the cohesion 
arising from the formation of liquid bridges between ad- 
jacent particles, the static as well as dynamical behav- 
ior of wet granular matter, such as the wet sand used for 
sculptures, differs from its dry counterpart dramatically 
[7, 8]. A former investigation on the pattern formation 
of vertically agitated wet granular layer reveals a pecu- 
liar period tripling spiral pattern, typically composed of 3 
arms that rotate continuously [1]. The preferred number 
of arms arise from period tripling of the agitated granular 
layer, which breaks the symmetry and drives the rotation 
of spiral arms. In the current work, the rotation frequency 
of the spiral arms is characterized quantitatively with two 
methods: One based on the horizontal motion of the arms 
from the top view snapshots, and the other one based on 



the fluctuations of the surface profile from the laser pro- 
filometry method. 

Figure 1(a) is a sketch of the main parts of the exper- 
imental setup. Cleaned glass beads (SiLiBeads S) with a 
diameter of d = 0.78 mm and 10% polydispersity, after 
proper mixing with purified water (Laborstar TWF-DI), 
are used as the granular sample. The amount of wetting 
is characterized by the liquid content W = V w /V g , where 
V w is the volume of the water and V g is that of the glass 
beads. It is normally kept within a few percent so that 
the cohesion arises mainly from the formation of cap- 
illary bridges between adjacent particles (as illustrated 
in the dash box). A certain amount of the sample (mass 
m = 1 13 g, corresponding roughly to 4-5 layers) is filled 
into a cylindrical polycarbonate container with an inner 
radius R = 8 cm and a height of H = 1 .06 cm. Sinusoidal 
vibrations generated by an electromagnetic shaker (Tira 
TV50350) are used to agitate the sample against grav- 
ity. The vibration frequency / and the nondimensional 
acceleration T = 47l 2 f 2 A/g, with vibration amplitude A 
and gravitational acceleration g, are the two control pa- 
rameters. The mobility of the patterns are captured by a 
high speed camera (IDT MotionScope M3) mounted on 
top of the container. The camera is synchronized with 
the shaker by a multi-pulse generator, in order to capture 
images at fixed phases of each vibration cycle. A more 
detailed description of the experimental setup could be 
found elsewhere [1]. 

Two methods are used to capture the mobility of the 
spiral arms. The first one uses the snapshots of the pattern 
obtained with a low angle illumination, which provides 
the information on the horizontal motion. The second 
one uses the the laser profilometry method [9], which 
provides information on the surface height fluctuations. 
As sketched in Fig. 1(b), the laser sheet illuminates a line 
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FIGURE 1. (a) A sketch of the experimental setup to capture 
the spiral pattern in a thin layer of wet granular matter under 
vertically sinusoidal agitations. The zoom-in view shown in 
the dotted box illustrates a capillary bridge formed between 
two adjacent particles, (b) A top view image captured with a 
combined illumination of a laser sheet and a LED ring, (j) is 
the angle between the laser sheet (illustrated as half transparent 
red) and the horizontal plane. Yellow dash curves highlight the 
spiral arms. 



of the sample surface with a fixed angle (j) « 37°, which 
transfers the height profile of the illuminated line into the 
images plane. 

From the top view images, the rotation frequency of 
the spiral arms can be obtained from the gray value 
fluctuations g(t) of regions with fixed distance to the core 
of the spirals, based on the fact that the spiral arms in the 
snapshots are brighter than the other regions, as Fig. 2(a) 
shows. The gray value difference is due to the fact that 
spiral arms correspond to kinks separating regions with 
different heights [1], and more light from the low angle 
illumination will be reflected to the camera in the kink 
(spiral arm) regions. 

The essential part of the algorithm is to choose the re- 
gion of interest for the analysis of the brightness fluctu- 
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FIGURE 2. Procedure to obtain rotation frequency from top 
view snapshots, (a) A typical image of a three armed rotating 
spiral pattern (averaged over 3 consequent snapshots to en- 
hance the contrast). The yellow cross is the current position of 
the spiral core determined by the image processing procedure, 
and the yellow curve corresponds to the trajectory of the spiral 
core in the past 12.5 seconds (1000 vibration cycles), (b) A 
close view of the spiral core region. The regions used to obtain 
gray value fluctuations marked with squares have equal dis- 
tance r to the spiral core, (c) Averaged power spectrum |g(v)| 2 
of the gray value fluctuations g(t) obtained within 6 regions. 
Inset is the fluctuation with time (blue dots) and a fitting with a 
harmonic function (red curve). 



ations properly. Especially the regions close to the spiral 
cores should be avoided, since those regions correspond 
to the common part of spiral arms, where the brightness 
fluctuations are much smaller than those further away 
from the core. Therefore, our first step is to locate the 
spiral cores. Based on the fact that the spiral core region 
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FIGURE 3. Procedure to obtain rotation frequency from 
laser profilometry. (a) A snapshot captured with laser sheet 
illumination, from which the surface profile (curve with half 
transparent red) is obtained by image processing, (b) Time- 
space plot of the surface profile captured with a frame rate 
f/3. h denotes the height of the profile, (c) Power spectrum of 
the height fluctuations |/z(v)| 2 , averaged over 6 positions with 
equal distances. The inset shows the height fluctuation with 
time at one position. 



keeps its brightness during the rotation of spiral arms, 
a minimization of the gray value over continuous snap- 
shots effectively enhances the contrast between the spiral 
core and the surrounding regions. The number of snap- 
shots is chosen to cover 1 /3 of the rotation period. After 
proper thresholding, we are able to estimate the location 
of the spiral core from the weighted average of the core 
region. 

From the trajectory of the spiral cores determined 
(see the superimposed curve in Fig. 2(a)), no clear re- 
lation to the rotation direction of the spiral arms could be 
found and the core wander around randomly. As shown 
in Fig. 2(b), 6 square boxes with a certain radial dis- 
tance r to the spiral core, and equally distributed in the 
angular direction, are chosen as the regions of interest. 
The spatially averaged brightness within each box as a 
function of time is recorded individually (see the inset 
of Fig. 2(c) for an example). As shown in Fig. 2(c), 
the averaged power spectrum of the brightness over all 
the 6 boxes indicates clearly a periodicity of the gray 



value fluctuations. From the peak value (|g(v p o)| 2 at 
frequency v p o) and the largest nearest neighbor value 
(|g(v p i)| 2 at frequency v p i), we interpolate the rotation 
frequency with a weighted average v r = (|g(v p o)| • v p o + 
l*(Vpi)| • v p i)/3(|g(Vpo)| + |s(Vpi)|), where the factor 3 
arises from the fact that all the three spiral arms will gen- 
erate brightness fluctuations. Alternatively, the frequency 
of g(t) is also determined directly by a sinusoidal fit of 
the data in the time domain (red curve in the inset of 
Fig. 2(c)), which agrees with the v r from frequency do- 
main within an error of 1%. 

With the laser sheet illumination, the rotation of spi- 
ral arm is represented by the kink propagation along the 
illuminated line. Therefore the rotation frequency of the 
spiral arms corresponds to the repetition rate of the kink 
propagation. As shown in Fig. 3(b), the time- space plot 
of the height profile clearly demonstrates this periodic- 
ity. The height profile of every third vibration cycle is 
recorded, in order to avoid the influence from the fre- 
quency of the surface height fluctuations (f/3) owing to 
period tripling. At a certain position, the plateau with 
maximal height (red color), which represents a region 
bounded by two spiral arms, repeatedly appears. The in- 
set of Fig. 3(c) shows one example of the height fluc- 
tuations with time h(t). Following the above analysis 
on the g(t), the rotation frequency is obtained from the 
power spectrum of the height fluctuations average over 
6 equally distributed positions (see Fig. 3(c)). Note that 
the height fluctuation is closer to a Heaviside step func- 
tion than to a sinusoidal fluctuation, thus fitting directly 
on h(t) is not appropriate. The smaller peak at about 2v r 
arises from the height fluctuation induced by the third 
spiral arm. 

Utilizing both methods, we are now ready to charac- 
terize the dynamics of rotating spirals. Qualitative obser- 
vations reveal that the emerging spiral arms as the thresh- 
old acceleration is reached will have a finite rotation fre- 
quency. Our main task is to characterize the dependence 
of the rotation frequency on various control parameters. 
Here we focus only on the case of single spiral core and 
leave the case of multiple spirals for further investiga- 
tions, because the interactions between spiral arms linked 
to different cores may lead to a variation of the rotation 
frequency. 

Figure 4 shows the rotation frequency v r as a func- 
tion of r for two vibration frequencies obtained from 
the top view snapshots taken at different time: shortly 
after changing control parameter T or with a delay of 
At & 300 s. For both vibration frequency / and T, v r cap- 
tured at different time agree with each other within the 
errorbar, suggesting that the rotation speed of the arms is 
time independent for a constant driving condition. As T 
is beyond the threshold value T t h, the rotation frequency 
increases monotonically with T within a certain range. 
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FIGURE 4. Rotation frequency of the spiral arms v r as a 
function of the peak vibration acceleration T for two vibration 
frequencies (J) with liquid content W = 1.5%. The threshold 
for spiral pattern to emerge (T t h) is 22.5 ± 1 .3 for / = 70 Hz and 
30.5 ± 1 .5 for / = 80 Hz. The rotation frequencies are obtained 
from the top view snapshots of the pattern, which are captured 
either shortly after varying T (with a delay time At < 30 s) or 
after a delay of At ^ 300 s. 
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FIGURE 5. A comparison between the rotation frequency v r 
obtained from top view snapshots and that from the laser pro- 
filometry method for / = 80 Hz. Inset shows the dependence of 
v r on the deviation of T from the threshold value T± obtained 
from top view snapshots for both / = 70 Hz and 80 Hz. Solid 
lines are linear fits to the data. 



Beyond the upper limit, the probability to have a single 
spiral core is too low to obtain reliable data points. 

Since the rotation speed of spiral arms is stable with 
time, a comparison of the results obtained with both 
methods described above is possible. Figure 5 shows 
such a comparison for / = 80 Hz and various T, ob- 
tained by collecting both g(t) and h(t) information con- 
sequently for each run of experiment. The error bars cor- 
respond to the variation of v r for different runs of the 
experiment. The good agreements between the results 
from both methods indicate that the dependence of v r 
on r obtained does not rely on the method chosen, thus 



either method can act as a counterproof of the other. 
Technically speaking, the analysis of the top view snap- 
shots, although it provides more straightforward infor- 
mation, will cost much more computing power compared 
with the analysis based on laser profilometry method, be- 
cause the latter contains only one dimensional informa- 
tion. Therefore the laser profilometry method will be the 
better candidate for the future analysis on a wide range 
of control parameters. 

In the inset of Fig. 5, the initial rotation frequency at 
the threshold acceleration v c (r t h) and the slope of v r 
with r for both frequencies are compared. The depen- 
dence of v r on r obtained so far suggests a linear depen- 
dence. Linear fits of the dependence on T — T t h yields a 
common threshold frequency v c = 0.45 Hz, and a larger 
slope of Av r /Ar = 0.11 Hz for / = 80 Hz, compared to 
the slope 0.05 Hz for 70 Hz. 

In conclusion, the rotation dynamics of spiral pattern 
in vertically agitated wet granular layer are characterized 
with two methods: one based on the horizontal motion of 
the spiral arms and the other one based on the vertical 
height fluctuation of a single illuminated line. For the 
case of single spiral core, both methods reveal that the 
rotation frequency of spiral arms grows linearly with the 
peak vibration acceleration and there exists a finite initial 
speed at the threshold acceleration where the pattern 
emerges. 

Detailed analysis on a broader range of driving param- 
eters, and for the case of multiple spiral cores will be a 
focus of further investigations. Since the time scale of 
the rotation arises presumably from the internal granu- 
lar temperature gradient [1], the analysis will help us to 
gain insight into the dynamics of fluidized wet granular 
matter. 
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